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The  kinetics  of  the  oxygen  reduction  reaction  (ORR)  occurring  on  nitrogen-doped  carbon  nanotubes 
(NCNT)  and  carbon  supported  platinum  (Pt/C)  catalysts  in  alkaline  conditions  is  analyzed  through  rotat¬ 
ing  ring  disk  electrode  (RRDE)  measurements  and  modeled  based  on  typical  mechanisms  reported  in  the 
literature.  The  purely  chemical  intermediate  reactions  are  found  to  be  insensitive  in  the  modeling  of  this 
mechanism,  whereby  the  models  can  be  simplified  to  the  occurrence  of  parallel  two  and  four  electro¬ 
chemical  electron  pathways.  Rotating  ring  disk  measurements  are  shown  to  be  an  alternative  technique 
to  split  the  total  current  density  into  partial  currents  that  could  be  associated  with  the  respective  reaction 
pathways.  Experimental  comparisons  of  the  NCNT  and  Pt/C  catalysts  revealed  that  the  ORR  kinetics  on 
these  materials  proceeds  similarly  throughout  two  and  four  electron  pathways.  The  four-electron  path¬ 
way  is  found  to  be  the  most  dominant  in  both  catalysts.  Kinetic  improvements  for  the  NCNT  catalysts  are 
observed,  positioning  them  as  promising  materials  for  the  ORR  catalysis  in  fuel  cells. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  the  demand  for  new  materials  for  low  tem¬ 
perature  fuel  cells  applications  has  increased  considerably.  This 
perpetuated  interest  stems  primarily  from  the  limited  performance 
and  high  system  cost  of  this  technology.  It  is  of  prime  concern  to 
develop  electrocatalyst  materials  with  improved  oxygen  reduc¬ 
tion  reaction  (ORR)  kinetics,  while  simultaneously  reducing  the 
excessive  cost  of  the  currently  utilized  carbon  supported  platinum 
(Pt/C)  metal  catalysts.  Pt/C  has  been  demonstrated  as  an  excellent 
material  for  ORR  catalysis;  however,  the  limited  natural  abundance 
and  steadily  increasing  cost  renders  its  long  term  applicability 
unpractical  for  commercial  purposes.  To  this  concern,  non-noble 
electrocatalysts  have  been  developed  as  alternatives  to  Pt-based 
catalysts  either  in  acidic  [1]  or  alkaline  conditions  [2-4].  In  alka¬ 
line  conditions,  several  different  types  of  non-noble  catalyst  have 
been  explored,  including  transition  metal  oxides  [3]  or  alloys  [4]; 
however,  their  current  state  of  development  is  not  sufficient  for 
application  in  new  age  fuel  cells.  Recently,  greatly  enhanced  ORR 
activity  of  CNT  has  been  attained  by  nitrogen  doping,  which  leads  to 
the  formation  of  nitrogen-doped  carbon  nanotubes  (NCNT).  NCNT 
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catalysts  with  better  electrocatalytic  activity  in  alkaline  conditions 
and  long-term  stability  compared  to  Pt/C  have  been  synthesized 
[5],  with  numerous  studies  elucidating  the  nature  of  performance 
improvements  and  essential  properties  of  these  materials  [6-8].  For 
NCNT,  nitrogen  on  the  surface  of  the  CNT  is  an  essential  require¬ 
ment  to  obtain  catalytically  active  sites. 

To  date,  comprehensive  modeling  to  account  for  the  ORR  mech¬ 
anism  occurring  on  these  catalysts  has  not  been  addressed.  A 
deeper  understanding  of  the  ORR  mechanism  occurring  on  NCNT 
would  allow  further  improvement  of  the  ORR  kinetics  occurring 
on  these  materials.  Therefore,  this  study  analyzes  the  ORR  mech¬ 
anism  on  highly  active  NCNT  non-noble  electrocatalysts  formed 
using  ethylenediamine  (EDA-NCNT)  as  a  nitrogen  and  carbon  pre¬ 
cursor  [1,6].  The  high  ORR  activity  observed  for  EDA-NCNT  and 
the  additional  benefits  of  using  CNT  based  materials  as  electrode 
catalysts  due  to  their  superior  mechanical  properties  [9],  thermal 
properties  [10],  and  low  cost,  positions  EDA-NCNT  as  a  potential 
replacement  material  to  Pt-based  catalysts  in  fuel  cell  cathodes. 
Experimental  measurements  using  a  rotating  ring  disk  electrode 
setup  are  used  to  evaluate  the  electrocatalytic  activity  of  these 
materials  along  with  commercial  Pt/C  for  comparison.  Three  differ¬ 
ent  physicochemical  models  accounting  for  the  diffusion,  as  well  as 
the  convection  of  adsorbed  and  dissolved  species  in  the  boundary 
layer  adjacent  to  the  working  electrode  are  developed  to  accurately 
describe  the  experimental  results.  These  models  are  based  on  typ¬ 
ical  reaction  mechanisms  reported  for  ORR  in  the  literature  under 
alkaline  conditions.  The  models  are  statistically  fit  to  the  obtained 
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experimental  data  in  order  to  estimate  the  kinetic  parameters  asso¬ 
ciated  with  each  model.  A  sensitivity  analysis  is  performed  for  the 
RRDE  measurements  to  discriminate  the  significance  of  the  reaction 
mechanisms  in  the  models,  and  to  elucidate  the  most  appropriate 
mechanism  for  NCNT  catalysis. 

2.  Experimental 

2.1.  NCNT  synthesis 

NCNT  synthesis  was  carried  out  by  a  simplistic  chemical  vapor 
deposition  (CVD)  growth  method.  Ferrocene  (98%,  Aldrich)  as  a 
growth  catalyst  was  dissolved  (2.5  wt.%)  in  ethylenediamine  (98% 
EMD  Chemicals)  precursor.  Small  quartz  tubes  were  utilized  as  the 
growth  substrate  and  were  placed  in  the  center  of  a  larger  alumina 
tube  running  through  the  center  of  a  furnace.  Synthesis  was  carried 
out  by  heating  the  tube  furnace  to  800  °C  and  injecting  the  precur¬ 
sor  solution  under  argon  protection.  Growth  was  carried  out  for  1  h, 
and  the  final  NCNT  products  were  obtained  from  the  inside  of  the 
quartz  tube  substrates.  NCNT  were  washed  in  0.5  M  sulfuric  acid 
for  5  h  and  dried  in  an  oven  overnight. 

2.2.  Oxygen  reduction  reaction  activity  evaluation 

The  electrocatalytic  activity  towards  ORR  was  evaluated  using 
RRDE  voltammetry.  Testing  was  carried  out  using  a  bipotentiostat 
(Pine  Instrument  Co.,  AFCBP-1 )  connected  to  a  rotation  speed  con¬ 
troller  (Pine  Instrument  Co.,  AFMSRCE).  Experiments  were  carried 
out  in  an  alkaline  electrolyte  using  100  mol  m-3  KOH.  A  Pt  wire 
counter  electrode  and  a  double  junction  Ag/AgCl  reference  elec¬ 
trode  were  utilized  for  all  testing,  and  the  potentials  referred  to  in 
this  work  are  on  this  scale.  A  glassy  carbon  electrode  surrounded 
by  a  Pt  ring  was  used  as  the  working  electrode.  20  p,L  of  a  cat¬ 
alyst  ink  solution  (4  mg  of  catalyst  dispersed  in  2  mL  of  0.2  wt.% 
Nation  solution)  was  directly  applied  onto  the  glassy  carbon  disk 
electrode.  The  working  electrode  potential  was  swept  from  0.2  to 
-1 .2  V  vs.  Ag/AgCl  (figures  only  to  -0.5  V  vs.  Ag/AgCl)  at  a  scan  rate 
of  0.01  Vs-1.  The  ring  potential  was  maintained  constant  at  0.5  V 
vs.  Ag/AgCl.  Before  testing,  the  electrolyte  was  saturated  with  nitro¬ 
gen  gas  and  a  potential  sweep  was  performed  in  order  to  obtain  a 
background  reading,  which  was  removed  from  the  subsequently 
measured  currents  obtained  during  ORR  testing.  After  the  back¬ 
ground  scan,  oxygen  gas  was  introduced  to  saturate  the  solution. 
The  experiments  were  performed  at  the  same  scan  rate  with  differ¬ 
ent  working  electrode  rotation  speeds.  A  saturated  02  environment 
was  kept  during  the  experiments,  which  ensured  the  collection  of 
reproducible  voltammograms. 

3.  Modeling 

The  models  used  in  the  present  study  consider  the  kinetics 
of  typical  ORR  processes  occurring  on  NCNT  and  Pt/C  catalysts. 
On  these  surfaces,  the  ORR  mechanism  is  complicated  and  could 
involve  several  different  intermediates.  The  identity  of  interme¬ 
diate  compounds  depends  primarily  on  the  composition  of  the 
catalyst  and  electrolyte.  In  this  section,  three  physicochemical 
models  are  developed  considering  the  kinetics  of  ORR  mechanisms 
that  have  been  previously  proposed  [11-16].  In  the  next  section, 
the  models  are  applied  to  analyze  the  RRDE  experimental  data 
and  evaluate  their  validity.  The  first  proposed  ORR  mechanism  is 
based  on  the  study  reported  by  Xu  et  al.  [14]  and  Appel  et  al.  [15]. 
In  this  reaction  mechanism,  elementary  steps  are  considered  for 
the  reduction  of  oxygen  (02)  and  hydroperoxide  (H02_)  species.  In 
the  last  reaction  of  the  proposed  mechanism,  the  H02_  generated 


by  two  different  pathways  reacts  with  water  molecules  to  form 
hydroxyl  ions  via  a  two-electron  transfer  reaction. 


O2  +  e  02ads  (1) 

02ads“  +  H20  H02+OH-  (2) 

H02  +  02-^  H02-  +  02  (3a) 

or 

H02  +  e_  H02-  (3b) 

H02-  +  H20  +  2e”  -*  30H-  (4) 


In  the  above  equations,  the  subindex  ads  stands  for  an  adsorbed 
state  of  the  species,  whereas  soluble  species  are  specified  without 
a  subindex.  The  second  and  third  models  investigated  are  based  on 
the  typical  ORR  mechanisms  involving  the  straightforward  two- 
and  four-electron  reduction  pathways.  The  kinetics  of  model  2 
only  considers  a  two-electron  pathway  where  molecular  oxygen  is 
reduced  to  hydrogen  peroxide.  This  intermediate  is  further  reduced 
to  form  OH-  [11-13,16]. 

02  +  H20  +  2e“  -*  H02-  +  OH-  (5) 

H02-  +  H20  +  2e"  -*  30H-  (6) 

The  third  model  entails  the  former  two-electron  pathways  (reac¬ 
tions  (5)  and  (6)),  in  addition  to  a  parallel  direct  four-electron 
pathway. 

02  +  2H20  +  4e"  -*  40H-  (7) 

More  detailed  ORR  mechanisms  have  been  reported  in  the  liter¬ 
ature  based  on  several  different  materials,  including  NCNT  catalysts 
[13,14,17,18].  Although  some  of  these  mechanisms  can  be  very 
complex  [17,18],  the  characterization  techniques  used  in  this  work 
(e.g.  voltammetry)  are  not  sensitive  to  most  of  these  reaction  steps. 
Possible  causes  of  the  lack  of  sensitivity  could  be  the  short  time 
domains  and  life-time  of  the  species  generated  on  the  surface  of 
the  catalysts.  Specifically,  their  generation-consumption  rates  are 
extremely  fast  compared  to  the  detection  limit  of  the  technique. 
Therefore,  only  those  reaction  mechanisms  where  the  detection 
limit  is  not  entirely  restricted  have  been  considered. 

3.1.  Electrode  reaction  mechanism  and  kinetics 

The  partial  current  density  for  electrochemical  reactions  (1), 
(3b)  and  (4)-(7)  considering  Langmuir  adsorption  behavior,  and 
without  assuming  rate-controlling  steps  are  described  below.  The 
backward  direction  of  each  step  is  neglected  on  these  expressions. 


Reaction  (1): 

J  =  ki  exP  (-“1  w)  (1  “  0°2,^~K2  (8) 

Reaction  (3b): 

Jjr  =  k3b  exp  (-a3b  ^  )  c*  o2  (9) 

Reaction  (4): 

^=k4exp(-2a45S)CV  (10) 

Reaction  (5): 

^=k5exp(-2a5^)cS2  (”) 

Reaction  (6): 

§=ksexp(-2^w)CH02-  (12) 
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Reaction  (7): 

ZF=k7exp(-4a7w)Co2  03) 

where  Cq2  ’  Cho2  and  Cho2_  rePresent  t^e  surface  concentrations  of 
02,  H02  and  H02_,  respectively,  while  0q2  is  the  fraction  of  the 
electrode  area  covered  by  02  ads-.  Moreover,  kx  and  ax  are  the  rate 
constant  and  charge  transfer  coefficient,  respectively  for  reaction  x 
and  F  is  the  electrode  potential  corrected  for  the  electrolyte  ohmic 
resistance  Rs ,  i.e., 

E'  =  E-jRs  (14) 


Table  1 

Constants  obtained  from  the  fitting  procedure  of  model  1. 

Parameter 

Value 

Do2 

3  x  10-9  m2  s_1 

Dho2 

2.45  x  10_1°  m2  s-1 

dho2~ 

2.65  x  10_1°  m2  s-1 

V 

1  x  10-6  m2  s_1  a 

Eho2~ 

1.0  x  10-5  m2  s_1 

a  Ref.  [23]. 


and  convective  transport  [19].  In  addition,  02  molecules  are 
considered  neutral. 


In  voltammetry,  the  potential  E  is  varied  linearly  with  time  at  a 
known  rate,  i.e., 

E  =  ^ini  —  Sratef  (15) 

where  Eini  is  the  starting  potential  (V)  and  srate  is  the  scan  rate 
(Vs-1). 

The  faradaic  current  density  for  model  1  is  given  by  the  following 
equation: 

JMI  =j\  +J*3f)  +J*4  (15) 

Model  2  involves  reactions  (5)  and  (6).  Thus,  the  faradaic  current 
density  for  this  model  yields: 

Jm2  —  J5  T J6  (17) 

Additionally  to  reactions  (5)  and  (6),  model  3  entails  a  direct  four 
electron  pathway  (reaction  (7)).  The  faradaic  current  density  of  this 
model  can  be  defined  as: 

Jm3  =h  +J6  +  J7  (18) 

The  total  current  density  measured  during  the  potential  scan  is  only 
composed  of  faradaic  ( jMy )  since  the  capacitive  component  has  been 
removed  from  the  experimental  measurements  (e.g.  background 
reading): 

jj  =  )My  (19) 


The  solution  is  separated  into  a  boundary  layer  region  of  thick¬ 
ness  S  and  a  well-mixed  bulk  region  where  concentrations  are 
uniform.  The  bulk  region  is  considered  to  be  fully  established  at 
a  distance  of  3 8  from  the  electrode  [20,21  ]: 

8  =  \.6\d]/3£2-V2vV6  (21) 

where  D;  is  the  diffusion  coefficient  of  any  species  involved  in 
a  transport-controlled  reaction,  Q  is  the  rotational  speed  of  the 
working  electrode  and  v  is  the  kinematic  viscosity  of  the  solution 
(see  Table  1  for  numerical  values  used  in  this  study). 

The  one-dimensional  diffusional  and  convective  flux  N;(y,  t)  of 
the  dissolved  species  i  (e.g.  02,  H02_)  can  be  expressed  as: 

0C* 

N^-Di^Z+VyQ  (22) 

where 

vy  =  -0.51023 Q3l2v~'[l2y2  (23) 

The  transient  mass  balances  for  02  and  H02_  species  within  the 
boundary  layer  region  0  <y  <  38  are: 

|=-VNi  (24) 

The  concentration  of  each  species  aty  =  38  is  at  its  bulk  level  C?, 
i.e., 


where  y  is  the  number  of  the  model. 

The  material  balance  for  the  sites  on  the  catalyst  surface  occu¬ 
pied  by  adsorbed  species  (e.g.  02  ads_)  is  described  below. 

d0n  -  /  FF'\ 

r  _ — =/<iexp  (  -cii  4=  ]  (1  -6U  -)Csn  -k20n  - 

°2,ads  dt  RTJ  °2,ads  J  °2  Z  °2,ads 

(20) 

where  Fq2  ad  -  represents  its  adsorption  density  required  to  com¬ 
pletely  fill  a  monolayer. 

3.2.  Transport  model  development 


Q(38,  t )  =  Cf  (25) 

The  boundary  conditions  at  the  electrode  surface  y  =  0  are  given 
below  depending  on  the  model: 


Model  1: 


D 


9Cq2 
°2  dy 


y= o 


h 

F 


^3a^HQ2 


^H02 


9CHo2 

dy 


y= o 


-  k20o 


2,  ads 


+  ^3a^HQ2 


(26) 

(27) 


The  following  assumptions  are  made  to  develop  the  transient 
kinetic-transport  equations  applicable  to  a  rotating  ring  disk  work¬ 
ing  electrode: 


D] 


ho2- 


dy 


y= o 


Model  2: 


—  ^3a^H02 


hb  U 
F  2F 


(a)  The  current  distribution  is  uniform  over  the  RRDE;  only  trans¬ 
port  in  the  y-direction  normal  to  the  electrode  surface  is 
considered. 

(b)  The  transport  and  kinetic  parameters  are  uniform  throughout 
the  system. 

(c)  The  solution  is  isothermal  at  298  K. 

(d)  Ideally  dilute  behavior  is  assumed  in  the  solution  (i.e.  activity 
coefficients  are  equal  to  1.0). 

(e)  Since  the  solution  contains  a  supporting  electrolyte,  migration 
of  electroactive  species  is  neglected  in  comparison  to  diffusional 


D 


dCo2 
°2  dy 


y= o 


h_ 

2  F 


D, 


ho2- 


^CH02- 

dy 


y= o 


Model  3: 


h_ 

2  F  2  F 


D0? 


dQ 


02 

dy 


y= o 


=  h_  h_ 

~  2F  4F 


(28) 


(29) 

(30) 


(31) 
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D, 


dC 


ho2- 


H02- 


dy 


y= o 


J-5  J6 

2F  2F 


(32) 


The  model  is  completed  by  specifying  that  the  concentration  of 
each  species  is  initially  uniform  everywhere  at  its  bulk  level  CP, 
i.e., 


c,(y,  0)  =  Cf  (33) 

while  the  surface  coverage  of  02iads_  is  zero  on  the  surface  of  the 
catalyst,  i.e., 

3.3.  Numerical  solution  of  model  equations  and  parameter 
estimation 


60 


E-jRs  /  V  vs.  Ag/AgCI 


Numerical  values  for  the  various  model  parameters  are  esti¬ 
mated  by  statistically  fitting  the  model  to  the  experimental  data. 
This  involves  the  combination  of  standard  methods  for  parameter 
estimation  with  numerical  solution  of  the  model  equations  pre¬ 
sented  in  the  previous  sections.  The  system  of  PDEs  given  for  each 
model  has  been  solved  using  the  finite  element  method  in  the  COM- 
SOL  Multiphysics™  4.1  software  package.  A  second-order  Lagrange 
quadratic  polynomial  is  used  as  the  shape  function.  An  interactive- 
smart  mesh  is  used  to  discretize  the  domain. 

The  model  is  statistically  fit  to  the  experimental  RRDE  data 
to  estimate  model  parameters  by  minimizing  the  sum-of-squares 
error  between  the  model  predictions  and  data.  The  fitness  function 
(35)  based  on  simplex  search  and  genetic  algorithms  provided  by 
the  MATLAB  R201  la  toolbox  was  used  [22]: 

Fitness  function  =  (jy  -j^xp)2  (35) 

where  jyT  and  j^xp  are  the  model-predicted  and  experimental  total 
current  densities,  respectively,  obtained  at  each  electrode  potential 
considered  along  a  scan.  An  Intel  Core  i7  CPU  running  at  2.80  GHz 
with  a  RAM  memory  of  8.00  GB  was  used  to  carry  out  the  calcula¬ 
tions. 

4.  Model  fitting  and  discussion 

A  3-fold  methodology  was  adopted  in  this  work  in  order  to 
analyze  the  influence  of  different  reaction  steps  in  the  overall 
mechanism  throughout  sensitivity  analysis.  The  first  step  was  to  fit 
the  experimental  data  obtained  for  the  NCNT  catalysts,  and  deter¬ 
mine  the  kinetic  parameters  accounting  for  the  ORR  mechanism.  As 
mentioned  previously,  electrochemical  techniques  including  lin¬ 
ear  sweep  voltammetry  are  not  sensitive  enough  to  consider  some 
reactions  due  to  the  short  time  domains  and  life-times  of  several 
species  generated  on  the  surface  of  the  catalysts.  Reactions  with 
time  durations  longer  than  the  discrete  detection  limit  were  deter¬ 
mined  in  the  second  step.  Finally  the  ORR  reaction  mechanism  was 
simplified  by  considering  only  the  elementary  reactions  sensitive  to 
the  experimental  measurements.  This  procedure  is  detailed  below. 

4A.  Analysis  of  rotating  disk  electrode  measurements  using 
models  t  and  2 

Fig.  1  shows  the  RRDE  voltammetry  used  to  evaluate  the 
electroacatalytic  activity  of  the  ORR  on  NCNT  catalysts.  These 
voltammograms  were  obtained  at  different  working  electrode  rota¬ 
tion  speeds  ( 1 00, 400, 900  and  1 600  RPM)  in  1 00  mol  nrr3  KOH  with 
oxygen  saturation  maintained  at  1  mol  m-3  by  continuous  bubbling 
[24].  A  sweep  rate  of  0.01  Vs-1  was  used  to  scan  the  potential 
from  0.2  V  vs.  Ag/AgCI  in  the  negative  direction.  The  low  over¬ 
potential  region  of  these  voltammograms  is  characterized  by  an 


Fig.  1.  Model-fitted  (black  continuous  line)  and  experimental  (symbols)  linear 
voltammograms  measured  at  0.01  V  s-1  and  different  rotation  speeds  on  glassy 
carbon  electrodes  coated  with  nitrogen-doped  carbon  nanotubes.  A  solution  of 
100  mol  m-3  KOH  was  used  as  supporting  electrolyte.  Fits  correspond  to  model  1. 
Experimental  data  reprinted  with  permission  from  Ref.  [2].  Copyright  2009  Ameri¬ 
can  Chemical  Society. 

independence  of  the  current  density  on  the  rotation  imposed  to  the 
electrode,  remarking  the  kinetic  control  of  ORR.  A  gradual  increase 
of  the  current  density  and  dependence  of  rotation  is  observed  to 
begin  at  about  -0.1  V,  suggesting  the  onset  of  the  mixed  control. 
At  more  negative  potentials,  a  limiting  current  density  is  attained 
for  all  rotations,  denoting  the  mass-transport  control  of  the  ORR.  At 
-0.1  V,  the  total  number  of  electrons  calculated  for  the  NCNT  and 
Pt/C  catalysts  were  3.63  and  3.55,  respectively.  This  implies  that 
the  ORR  on  these  surfaces  entails  more  than  3  electrons,  indicating 
that  more  than  two  reactions  might  be  implicit. 

Continuous  lines  in  Fig.  1  indicate  the  computed  polarization 
curves  obtained  by  the  fit  process  at  different  rotation  rates.  A 
good  agreement  is  observed  between  the  experiments  and  the 
model  predictions.  Slight  deviations  are  obtained  for  the  mixed 
region,  which  could  be  associated  with  additional  contributions  not 
considered  by  the  model,  e.g.  specific  adsorption  of  ions  [16]  and 
heterogeneities  of  the  catalyst.  The  constants  and  kinetic  parame¬ 
ters  obtained  from  this  procedure  are  reported  in  Tables  1  and  2, 
respectively.  It  is  worth  mentioning  that  an  effort  was  carried  out 
to  get  these  constants  from  literature  with  limited  success.  There¬ 
fore,  they  were  incorporated  in  the  fitting  sub-routine  along  with 
the  kinetic  parameters.  It  was  found  that  the  diffusion  coefficients 
and  the  maximum  adsorption  density  (refer  to  Table  1)  are  in  the 
range  of  typical  diffusion  coefficients  reported  for  other  species 
with  similar  atomic  size  [20,21,23]. 

A  sensitivity  analysis  was  performed  to  determine  the  signif¬ 
icance  of  the  reactions  steps  involved  in  the  ORR  mechanism  of 
model  1.  The  fit  of  the  model  is  very  sensitive  to  some  parameters 
and  insensitive  to  others.  To  more  closely  assess  the  certainty  of 
these  estimates,  simulations  were  performed  by  decreasing  one 
parameter  by  10%  while  keeping  all  others  fixed  at  the  values 


Table  2 

Parameters  obtained  from  the  fit  of  model  1  to  the  experimental  data  obtained  on 
NCNT  catalysts. 


Parameter 

Value 

OL\ 

0.47 

a4 

0.42 

Oi  5 

0.21 

ki 

5.85  x  10-5  ms-1 

/<2 

10  m2  s_1 

k3 

8  x  10-11  ms-1 

k4 

4.91  x  10-5  ms-1 

k5 

1  x  10-5  ms-1 
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Fig.  2.  Computed  partial  and  total  current  densities  for  the  electrochemical  reac¬ 
tions  considered  in  model  1. 

obtained  from  the  initial  model  fitting.  The  sensitivity  of  the  elec¬ 
trode  response  to  changes  in  variables  or  parameters  associated 
with  a  particular  step  of  a  reaction  mechanism  also  reflects  the 
extent  to  which  the  overall  reaction  rate  is  controlled  by  that 
step.  The  electrode  response  is  sensitive  to  changes  in  variables  or 
parameters  associated  with  the  rate-controlling  step.  Conversely,  if 
a  step  does  not  control  the  rate,  then  the  electrode  response  will  be 
insensitive  to  changes  associated  with  it.  The  results  of  this  analy¬ 
sis  were  very  straightforward  -  only  those  reaction  steps  involving 
charge  transfer  are  sensitive  to  the  model,  whereas  the  chemical 
reaction  (2)  between  superoxide  (02ads_)  and  water,  as  well  as 
the  regeneration  of  oxygen  (3a)  has  fast  kinetic  rates  beyond  the 
detection  limit  of  electrochemical  measurements.  The  polarization 
curves  in  Fig.  2  reflect  this  conclusion  where  the  total  current  den¬ 
sity  can  be  split  into  three  electrochemical  reactions.  According  to 
this  figure,  two  electrons  are  transferred  in  the  first  two  electro¬ 
chemical  steps  (reactions  (1)  and  (3b)),  followed  by  the  reduction 
of  hydroperoxide  species  (H02_)  which  involves  the  transfer  of 
two  additional  electrons  (reaction  (4)).  Thus,  as  considered  in  other 
studies  [11-13],  the  reaction  mechanism  can  be  simplified  to  two 
electrochemical  steps  each  involving  two-electron  transfer  (i.e. 
model  2). 

Similar  fitting  procedures  were  applied  for  model  2.  Fitting 
model  2  to  the  experimental  data  recorded  for  NCNT  catalysts 
yielded  similar  results  compared  to  Fig.  1,  and  therefore,  they  are 
not  shown  herein.  The  kinetic  parameters  accounting  for  the  ORR 
mechanism  are  detailed  in  Table  3.  As  expected,  the  diffusion  coef¬ 
ficients  of  02  and  H02_  species  are  similar  to  those  found  by  model 
1  (refer  to  Table  1).  The  sensitivity  analysis  carried  out  with  the 
kinetic  parameters  of  model  2  revealed  that  reactions  (5)  and  (6) 
are  sensitive  to  the  model,  and  this  reaction  mechanism  accounts 
accurately  for  the  ORR. 

Further  insight  into  the  system  was  provided  through  the 
computation  of  non-measured  variables  (e.g.  surface  concentra¬ 
tions).  Fig.  3  shows  the  effects  of  two  different  rotations  (100 
and  1600RPM)  on  the  surface  concentrations  (Q)  of  02  and  H02_ 
species.  According  to  this  figure,  C£  starts  to  decay  at  more  positive 


Fig.  3.  Variation  of  surface  O2  and  HO2-  concentration  during  the  potential  scans 
operated  at  0.01  V  s_1  and  two  different  rotations  ( 1 00  and  1 600  RPM).  The  calcula¬ 
tions  were  carried  out  using  model  2  and  the  parameters  shown  in  Table  3. 

potentials  when  the  rotation  applied  to  the  electrode  is  lower. 
This  effect  is  not  surprising  considering  that  lower  rotations  lead 
to  limited  convective  inputs  of  02  to  the  surface  of  the  catalyst. 
Thus,  the  depletion  of  02  occurs  earlier  for  the  experiment  carried 
out  at  100  RPM.  This  can  be  observed  in  Fig.  3  at  potentials  of  ca. 
-0.25  V,  and  -0.35  V  for  this  rotation  value.  These  potential  values 
are  located  where  the  input  of  mass-transport  control  is  observed 
(see  Fig.  1).  The  surface  concentrations  of  H02_  species  are  also 
shown  in  Fig.  3.  These  concentration  profiles  closely  follow  the 
surface  concentrations  of  02.  As  02  reduction  commences,  C*  _ 

HU2 

is  produced  at  a  constant  rate  until  the  reduction  of  02  becomes 
mass-transport  controlled.  At  this  point,  a  decline  in  C*  _  occurs 

HU2 

until  it  approaches  to  zero  at  more  negative  potentials.  Once  the 
depletion  of  H02_  species  begins,  the  reduction  of  02  and  H02- 
proceed  at  different  rates  evidenced  by  the  lower  02  concentra¬ 
tion  profile  throughout  this  region.  This  phenomenon  suggests  that 
the  incomplete  reduction  of  oxygen  (reaction  (5))  proceeds  at  a 
faster  rate  compared  to  the  reduction  of  the  hydroperoxide  (reac¬ 
tion  (6)).  Therefore,  H02-  species  are  not  consumed  as  fast  as  they 
are  produced. 

4.2.  Hydroxide  selectivity  ORR  mechanism  on  NCNT  catalysts 

In  the  previous  section,  two  different  models  were  analyzed  to 
determine  the  sensitivity  of  the  reactions  in  the  ORR  mechanism. 
It  was  determined  that  the  linear  sweep  voltammetry  technique  is 
sensitive  enough  to  measure  accurately  the  two  electron  pathway 
proposed  in  the  literature  [10-13].  Flowever,  an  important  question 
arising  at  this  point  is  whether  parallel  steps  in  the  ORR  mechanism 
could  occur  and  contribute  to  oxygen  reduction. 

Hydroxide  selectivity  in  alkaline  electrolytes  is  an  indicator  of 
the  efficiency  of  the  catalysts,  whereby  it  provides  insight  regard¬ 
ing  the  reduction  mechanisms  the  molecular  oxygen  during  the 
electrochemical  reaction  [16].  Hydroxide  selectivity  is  anomalous 
to  H20  selectivity  in  acidic  electrolytes  and  can  be  calculated  using 
the  following  equation: 


Table  3 

Parameters  obtained  from  the  fit  of  model  2  to  the  experimental  data  obtained  on 
NCNT  catalysts. 


Parameter 

Value 

a6 

0.32 

Oil 

0.18 

k6 

4.55  x  10-5  ms-1 

k7 

3.34  x  10-9  ms-1 

5OH — 


f  Jd  -  Ur/N)\ 
\jp  +  Or/AO  J 


x  100 


(36) 


where  j’d,  Jr  and  N  are  the  disk  current  density,  ring  current  density 
and  collection  coefficient  number  (0.26).  As  proposed  by  Yea¬ 
ger  et  al.,  ORR  can  occur  through  a  four  electron  pathway  where 
molecular  oxygen  is  reduced  to  form  OH-  (reaction  (7))  or,  a  two 
electron  pathway  where  molecular  oxygen  is  reduced  to  hydrogen 
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Fig.  4.  Model-fitted  (black  continuous  line)  and  experimental  (symbols)  linear 
voltammograms  measured  at  0.01  V  s-1  and  different  rotation  speeds  on  glassy  car¬ 
bon  electrodes  coated  with  nitrogen-doped  carbon  nanotubes  (main  plot)  and  a  Pt 
ring  (inset).  A  solution  of  100  mol  m-3  KOH  was  used  as  supporting  electrolyte.  Fits 
correspond  to  model  3  (refer  to  Table  4). 

Experimental  data  reprinted  with  permission  from  Ref.  [2].  Copyright  2009  Ameri¬ 
can  Chemical  Society. 

peroxide  (reaction  (5)),  which  further  reduces  to  form  OH-  (reac¬ 
tion  (6))  [12,13]. 

To  assess  the  simultaneous  occurrence  of  these  parallel  reac¬ 
tions,  we  conducted  rotating  ring  disk  electrode  experiments  by 
applying  a  reductive  potential  in  the  disk  and  an  oxidative  potential 
in  the  ring.  The  intention  of  these  measurements  was  to  identify  and 
quantify  the  soluble  species  generated  from  ORR  on  the  disk  by  the 
anodic  ring  electrode.  The  oxidative  flux  of  species  on  the  ring  can 
be  associated  with  the  hydroperoxide  produced  from  ORR  (reaction 
(5))  that  is  not  further  reduced  on  the  surface  of  the  electrocatalyst 
(reaction  (6))  and  is  released  into  the  electrolyte.  Consequently, 
the  current  density  of  the  disk  involves  the  current  associated  with 
reactions  (5)-(7)  occurring  on  the  electrocatalyst  materials. 

Fig.  4  shows  two  different  sets  of  linear  sweep  voltammograms 
where  the  currents  of  the  disk  and  the  ring  have  been  separated 
into  contributions  of  reaction  (5)  (inset  of  the  plot)  and  reactions  (6) 
and  (7)  (main  plot).  The  voltammograms  were  obtained  at  different 
rotations  ( 1 00, 400, 900  and  1 600  RPM)  of  the  electrode  in  the  elec¬ 
trolyte  solution  containing  100  mol  m-3  KOH  and  1  molrn-3  02.  A 
sweep  rate  of  0.01  V  s-1  was  used  to  scan  the  potential  from  the 
OCP  in  the  negative  direction.  A  comparison  of  the  current  den¬ 
sities  displayed  in  the  main  plots  and  those  shown  in  the  inset 
demonstrates  that  oxygen  is  mainly  reduced  through  the  4  electron 
pathway.  Experimental  evidence  of  this  finding  has  been  previously 
reported  in  literature  [2,16].  The  effects  of  mass-transport  control 
are  also  observed  in  Fig.  4  for  both  reaction  pathways  (reactions 
(5)  and  (7)).  At  high  overpotential,  a  plateau  in  current  density  is 
attained  as  a  result  of  the  oxygen  depletion  on  the  surface  of  the 
disk  electrode.  This  effect  also  impacts  the  oxidative  current  of  the 
ring  since  the  oxidation  of  H02~  depends  on  its  diffusion  subse¬ 
quent  to  oxygen  reduction.  Fig.  4  demonstrates  the  fitting  of  the 
experimental  data  by  model  3  (continuous  lines).  In  order  to  per¬ 
form  this  task,  the  fitness  function  used  for  models  1  and  2  had  to 
be  modified  to  account  for  the  partial  current  densities  from  the 
parallel  reactions: 

Fitness  function  =  (j!™dcl3  +jfode'3  -jfp  -  j“p)2 

+  X)  (J™°del3  ~j6XP)2  07) 

The  model  predictions  shown  in  Fig.  4  are  good  considering 
that  a  single  set  of  parameters  account  for  different  experiments 


Table  4 

Parameters  obtained  from  the  fit  of  model  3  to  the  experimental  data  obtained  on 
NCNT  catalysts. 


Parameter 

Value 

Oi\ 

0.13 

0i4 

0.32 

Ois 

0.14 

/<! 

5  x  10-6  ms-1 

k2 

1  x  10-7  ms-1 

k  3 

2.31  x  10-7  ms-1 

obtained  under  different  conditions.  In  addition,  the  heterogeneity 
(i.e.  anisotropy)  of  these  materials  is  an  important  effect  not  consid¬ 
ered  by  the  model.  The  most  significant  deviation  of  the  model  is  the 
prediction  of  the  current  plateau  related  with  reaction  (6)  (refer  to 
the  inset  of  Fig.  4).  This  deviation  is  not  surprising  since  that  current 
density  is  calculated  indirectly  using  an  oxidation  current  rather 
than  a  reduction  current.  Therefore,  the  model  is  not  capable  of 
identifying  those  intermediate  contributions,  which  H02-  species 
overcome  before  being  oxidized  on  the  ring,  e.g.  diffusion  and  con¬ 
vection  from  disk  to  ring.  However,  these  inputs  are  well  described 
before  mass-transfer  control  is  reached.  The  kinetic  parameters 
obtained  from  the  fitting  procedure  are  reported  in  Table  4.  It  was 
found  that  the  diffusion  coefficients  (Dq2  and  DH o2 )  did  not  vary 
with  respect  to  those  shown  in  Table  1  for  models  1  and  2  (e.g. 
these  parameters  rely  on  the  mass-transport  control  region).  How¬ 
ever,  the  kinetic  parameters  obtained  with  models  1  and  2  present 
significant  changes,  whereby  rotating  ring  disk  measurements  are 
needed  to  describe  properly  the  ORR  mechanism  of  NCNT  catalysts. 

4.3.  Hydroxide  selectivity  ORR  mechanism  on  Pt-catalysts 

For  conventional  carbon  electrodes  in  acidic  and  alkaline  con¬ 
ditions,  it  has  been  reported  that  the  ORR  generally  follows  the 
two  electron  pathway  (e.g.  reactions  (5)  and  (6))  [20].  In  order  to 
explore  this  occurrence,  ORR  was  conducted  onto  commercial  Pt/C 
catalysts,  and  thus,  its  kinetics  was  compared  to  that  of  the  NCNT 
catalysts.  Fig.  5  shows  experiments  carried  out  on  commercial  Pt/C 
surfaces.  As  observed  from  the  comparison  between  Figs.  4  and  5, 
the  current  densities  associated  with  the  reduction  processes  are 
higher  on  the  NCNT  catalysts  than  on  Pt/C.  These  differences  are  not 
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Fig.  5.  Model-fitted  (black  continuous  line)  and  experimental  (symbols)  linear 
voltammograms  measured  at  0.01  V  s-1  and  different  rotation  speeds  on  commer¬ 
cial  Pt-catalysts  (main  plot)  and  a  Pt  ring  (inset).  A  solution  of  100  mol  m-3  KOH  was 
used  as  supporting  electrolyte.  Fits  correspond  to  model  3  (refer  to  Table  5). 
Experimental  data  reprinted  with  permission  from  Ref.  [2].  Copyright  2009  Ameri¬ 
can  Chemical  Society. 


J.  Vazquez- Arenas  et  al.  /  Journal  of  Power  Sources  205  (2012)215-221 


221 


Table  5 

Parameters  obtained  from  the  fit  of  model  3  to  the  experimental  data  obtained  on 
Pt-C  catalysts. 


Parameter 

Value 

OL\ 

0.13 

a4 

0.30 

Oi  5 

0.13 

/<! 

3.3  x  10-6  ms-1 

k2 

8.6  x  10-8  ms-1 

/<3 

3.9  x  10-7  ms-1 

significant  but  they  emphasize  the  potential  application  of  NCNT  to 
replace  Pt/C  catalysts  in  fuel  cells.  The  qualitative  similarities  pre¬ 
sented  between  these  catalysts  might  suggest  a  similar  ORR  kinetic 
on  these  materials.  The  existence  of  a  low  current  density  on  the 
ring  (see  inset  in  Fig.  5)  confirms  this,  and  justifies  the  use  of  model  3 
to  account  for  the  experimental  data  shown  in  Fig.  5.  In  other  words, 
the  ORR  mechanism  on  Pt/C  surfaces  also  proceeds  through  4  and  2 
electron  pathways  (reactions  (5)-(7)).  Table  5  reveals  that  there  are 
not  significant  variations  with  the  kinetic  parameters  found  for  the 
NCNT  catalysts  (Table  4).  Flowever,  as  observed  in  Figs.  4  and  5,  the 
current  density  associated  with  the  two  and  four-electron  path¬ 
ways  for  the  ORR  mechanism  are  slightly  higher  onto  the  NCNT 
catalysts  (see  Fig.  4)  with  a  direct  comparison  provided  in  a  pre¬ 
vious  report  [2].  Not  surprisingly,  this  higher  activity  is  related  to 
the  higher  kinetic  parameters  (rate  constants  and  transfer  coeffi¬ 
cients)  obtained  on  NCNT  (Table  4).  Thus,  on  this  material  the  rate 
constants  of  the  two-electron  pathway  are  higher,  as  well  as  the 
transfer  coefficient  for  the  four-electron  pathway. 

5.  Conclusions 

ORR  kinetics  occurring  on  NCNT  and  Pt/C  catalysts  were  ana¬ 
lyzed  through  typical  mechanisms  reported  in  literature.  It  was 
found  that  not  all  reactions  involved  in  the  reaction  mechanisms 
are  sensitive  to  the  current  density  response.  Particularly,  those 
steps  involving  chemical  reactions  were  found  to  be  insensitive  in 
the  model,  as  a  result  of  their  rapid  occurrence.  On  the  other  hand, 
linear  sweep  voltammetry  was  sensitive  to  two  and  four  electron 
pathways  for  ORR.  The  kinetic  parameters  obtained  from  fitting  the 
models  to  experimental  data  strongly  depend  on  the  reactions  pro¬ 
posed  in  the  mechanism.  Thus,  the  insensitivity  of  some  reaction 
steps  in  the  model  and  the  excess  of  degrees  of  freedom  between 
the  experiments  and  the  model  predictions  can  lead  to  misestimate 
the  kinetic  parameters.  Rotating  ring  disk  measurements  showed 
to  be  an  alternative  technique  to  split  the  total  current  density  into 
partial  currents  that  could  be  associated  with  the  respective  reac¬ 
tion  pathways.  Thus,  the  determination  of  the  kinetic  parameters 
of  the  ORR  mechanism  was  improved  by  reducing  the  number  of 
degrees  of  freedom.  The  sensitivity  analysis  conducted  to  evaluate 
the  validity  of  the  different  ORR  mechanisms  analyzed  in  this  work 
showed  that  the  chemical  reaction  between  superoxide  and  water, 
as  well  as  the  regeneration  of  oxygen  in  the  first  model  presented 


fast  kinetic  rates  beyond  the  detection  limit  of  electrochemical 
measurements.  Thus,  the  kinetics  of  the  ORR  reaction  mechanism 
was  simplified  to  two  (model  2)  and  three  (model  3)  electrochem¬ 
ical  steps.  Both  models  considered  a  two-electron  pathway  where 
molecular  oxygen  is  reduced  to  hydrogen  peroxide,  whereas  model 
3  additionally  entailed  a  parallel  direct  four-electron  pathway. 
Although  it  was  determined  that  the  ORR  mechanism  on  NCNT  and 
Pt-catalysts  was  sensitive  to  the  two-electron  pathway  considered 
in  model  2,  a  comparison  of  NCNT  and  Pt/C  catalysts  revealed  that 
the  ORR  kinetics  on  these  materials  proceeded  via  both  two  and 
four  electron  pathways  in  parallel  as  described  in  model  3.  The  four- 
electron  pathway  was  found  to  be  the  most  dominant  and  model  3 
the  most  suitable  ORR  mechanism  occurring  on  these  catalysts.  A 
slight  kinetic  improvement  onto  the  NCNT  catalysts  was  observed 
based  on  half-cell  electrochemical  testing  and  by  determination  of 
the  ORR  kinetic  parameters  on  these  materials,  rendering  NCNT  as 
promising  materials  for  the  ORR  catalysis  in  fuel  cell  technology. 
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